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Abstract 

We derive a high-resolution formula for the L 2 -quantization errors of Riemann-Liouville 
processes and the sharp Kolmogorov entropy asymptotics for related Sobolev balls. We describe 
a quantization procedure which leads to asymptotically optimal functional quantizers. Regular 
variation of the eigenvalues of the covariance operator plays a crucial role. 
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1 Introduction 

Functional quantization of stochastic processes can be seen as a discretization of the path-space 
of a process and the approximation (coding) of a process by finitely many deterministic functions 
from its path-space. In a Hilbert space setting this reads as follows. 

Let (H,< •, • >) be a separable Hilbert space with norm || • || and let X : (£l,A, IP) — ► H be 
a random vector taking its values in H with distribution Px- For n 6 N, the L 2 -quantization 
problem for X of level n (or of nat-level log n) consists in minimizing 



1/2 

a£a " " / " a£a 



Emin \\X — a\\ 2 ) = || min ||X — d|| ||z, 2 (P) 



over all subsets a C H with card(a) < n. Such a set a is called n-codebook or n-quantizer. The 
minimal rath quantization error of X is then denned by 

e n (X) := inf ((Emin \\X - a\\ 2 ) 1/2 : a C H, card(a) < n] . (1.1) 

^ ago J 

Under the integrability condition 

E||X|| 2 <oo (1.2) 

the quantity e n (X) is finite. 
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For a given n-quantizer a one defines an associated closest neighbour projection 

TT a := J2 al C a (a) 

and the induced a-quantization (Voronoi quantization) of X by 

X a :=7T a (X), (1.3) 
where {C a (a) : a S a} is a Voronoi partition induced by a, that is a Borel partition of H satisfying 

C a {a) C V a (a) := {x £ H : \\x - a\\ = min \\x - b\\} (1.4) 

for every a £ a. Then one easily checks that, for any random vector X : — ► a C H , 

E||X-X'|| 2 > E\\X -X a \\ 2 = E min||A-a|| 2 

so that finally 

e n (X) = inf{(E||A - X\\ 2 ) 1/2 : X = f(X),f : H -» # Borel measurable, (1.5) 
card(/(H)) < n} 

= inf {(E||X - Xf) l/2 : X : ft -> H random vector, card(X(^)) < n} . 

Observe that the Voronoi cells V a (a),a S a are closed and convex (where convexity is a charac- 
teristic feature of the underlying Hilbert structure). Note further that there are infinitely many 
a-quantizations of X which all produce the same quantization error and X a is P-a.s. uniquely 
defined if Px vanishes on hyper planes. 

A typical setting for functional quantization is H = L 2 ([0, l],dt) but is obviously not restricted 
to the Hilbert space setting. Functional quantization is the natural extension to stochastic processes 
of the so-called optimal vector quantization of random vectors in H = WL d which has been extensively 
investigated since the late 1940's in Signal processing and Information Theory (see [1J, |S]). For the 
mathematical aspects of vector quantization in M. d , one may consult [3], for algorithmic aspects see 
[15. and "non-classical" applications can be found in |14| . |16j . For a first promising application 
of functional quantization to the pricing of financial derivatives through numerical integration on 
path-spaces see [T7] . 

We address the issue of high-resolution quantization which concerns the performance of n- 
quantizers and the behaviour of e n (X) as n — > oo. The asymptotics of e n (X) for M rf -valued random 
vectors has been completely elucidated for non-singular distributions Px by the Zador Theorem 
(see (3) and for a class of self-similar (singular) distributions by (Hj- In infinite dimensions no such 
global results hold, even for Gaussian processes. 

It is convenient to use the symbols ~ and ~, where a n ~ b n means a n /b n — > 1 and a n ~ b n 
means limsupn^^ a n /b n < 1. A measurable function ip : (s,oo) — ► (0, oo) (s > 0) is said to be 
regularly varying at infinity with index b G R if, for every c > 0, 

,• ip(cx) b 
hm — -— = c . 

x^oo ip[x) 

Now let X be centered Gaussian. Denote by Kx C H the reproducing kernel Hilbert space 
(Cameron-Martin space) associated to the covariance operator 

C x :H^H,C x y:=E(<y,X>X) (1.6) 
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of X. Let Ai > A2 > • • • > be the ordered nonzero eigenvalues of C x and let {uj : j > 1} be the 
corresponding orthonormal basis of supp(Px) consisting of eigenvectors (Karhunen-Loeve basis). 

( d \ 

If d : = dimi^x < 00, then e n (X) = e n ® N(0, Xj) , the minimal nth L 2 -quantization error of 

V =1 / 

d 

N(0, Xj) with respect to the ^-norm on M. d , and thus we can read off the asymptotic behaviour 
3=1 

of e n (X) from the high-resolution formula 

d , \ l/2d / W _|_ O \ (rf+2)/4 

e n (®iV(0,A i ))~g(d)v^F(nJ =1 A J -) n" 1 ^ as n ^ 00 (1.7) 

i=i 

where q(d) E (0, 00) is a constant depending only on the dimension d (see [3])- Except in dimension 
d = 1 and d = 2, the true value of q(d) is unknown. However, one knows (see 5 ) that 

/ d \ 1/2 

g(d) ~ as d -> 00. (1.8) 

\27re/ 

Assume dimi^x = 00. Under regular behaviour of the eigenvalues the sharp asymptotics of 
e n (X) can be derived analogously to (|1.7j) . In view of (|1.8j) it is reasonable to expect that the 
limiting constants can be evaluated. The recent high-resolution formula is as follows. 

Theorem 1 Let X be a centered Gaussian. Assume Xj ~ <p(j) as j — > 00, where ip : 

(s, 00) — > (0,oo) is a decreasing, regularly varying function at infinity of index —b < —1 for some 
s > 0. Set, for every x > s, 

i/)(x) — 1 

X(f(x) 

Then 

( {b\ h ~ l b \ 1/2 
e n (X) ~ I I - ] tp{\ogn)~ 1 / 2 as n^oo. 



A high-resolution formula in case b = 1 is also available (see jllj). Note that the restriction 

00 

—b < —1 on the index of (p is natural since Yl -\? < 00 • The minimal L r -quantization errors of 

i=i 

X, < r < 00, are strongly equivalent to the L 2 -errors e n (X) (see j2]) and thus exhibit the same 
high-resolution behaviour. 

A related quantization problem is the Kolmogorov metric entropy problem for the closed unit 
ball 

U x := jx G K X : \\x\\ Kx < 1} = {x G supp(P x ) : £ 1^^! < i| (1.9) 

of (Strassen ball). Note that Ux is a compact subset of H. For n e N, the metric entropy 
problem for {Jx consists in minimizing 

max min llx — all = II min ||X — alNLoofp-) 

over all subsets a C H with card(a) < n, where X is any iif-valued random vector with supp(P^' ) = 
Ux- The nth entropy number is then defined by 

e n (Ux) '■= inf < max min llx — all : a C H , card(a) < n > . (1-10) 
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If d := dim if x < oo, then e n (Ux) = e n (£d), the nth entropy number of the ellipsoid 

£ d := < x G 




with respect to the ^-norm on M d . Thus we can read off the asymptotic behaviour of e n (Ux) from 
the formula 

e n (£) ~ pCdJCn^AjO^Cvol (S d (0, l)))V«* n -i/«i as n ^ oo (1.11) 

where the constant p(d) 6 (0, oo) is unknown for d > 3 and p(d) ~ q(d), d — ► oo (see |S])- 

If dimi^x = oo, the recent solution of the Kolmogorov metric entropy problem for Ux is as 
follows. 

Theorem 2 (11^ ) Assume the situation of Theorem 1. Then 

'b\ b ' 2 



e n (Ux) ~ l-j if(logn) 1/2 as 



n — > oo. 



This formula is still valid for 6 = 1 and, ignoring the probabilistic interpretation, also for b > 
(00 := 1) provided Xj — > as j — * oo. (see [7], [T2]). A different approach via the inverse of e n (Ux), 
the Kolmogorov e-entropy, is due to Donoho 3 . (However, his result does not provide the correct 
constant.) ^,From Theorems 1 and 2 we conclude that functional quantization and metric entropy 
are related by 

/2 loe n\ x l 2 

e n {X) ~ ( y^rj en(Ux) as n ^ oo. (1.12) 

The paper is organized as follows. In Section 2 we investigate Riemann-Liouville processes in 
H = L 2 ([0, 1], dt). For p £ (0, oo), the Riemann-Liouville process X p = (A" t p ) tg [ 0j i] on [0,1] is 
defined by 

/"* i 

X p t := (t- s) p -2dW s (1.13) 



o 

where W is a standard Brownian motion. We derive a high-resolution formula for X p and cor- 
respondingly, the precise entropy asymptotics for fractional Sobolev balls. As a consequence we 
obtain a new result for fractionally integrated Brownian motions. In Section 3 we describe a quan- 
tization procedure which furnishes asymptotically optimal quantizers in the situation of Theorem 
1. Here the Karhunen-Loeve expansion plays a crucial role. In Section 4 we discuss a dimension 
conjecture. 



2 Riemann-Liouville processes 

Let X p = {X% )t e [o,i] be the Riemann-Liouville process of index p £ (0, oo) as defined in (|1.13|) . Its 
covariance function is given by 

EX p Xf= / (t-r) p -2( s -r) p -2dr. (2.1) 
Jo 

Using p A ^-Holder continuity of the application t \—* X p from [0,1] into L 2 (P) and the Kolmorogov 
criterion one checks that X p has a pathwise continuous modification so that we may assume without 
loss of generality that X p is pathwise continuous. In particular, X p can be seen as a centered 
Gaussian random vector with values in 

H = L 2 ([0,l],dt). 
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The following high-resolution formula relies on a theorem by Vu and Gorenflo |18j on singular 
values of Riemann-Liouville integral operators 

R f> 9® = tTo\ A* " sf'^&ds, 13 E (0, oo). (2.2) 



r(/3) Jo 

Theorem 3 For every p £ (0, oo), 

e n (X") ~ 7r-^\p + l/ 2 )P(^±l)V2 r( + 1/2)(1 

2p 

Proof. For /5 > 1/2, the Riemann-Liouville fractional integral operator Rg is a bounded operator 
from L 2 ([0, l],dt) into £ 2 ([0, l],dt). The covariance operator 

C p : L 2 ([0,l],df) -> L 2 ([0,l],di) 

of X p is given by the Fredholm transformation 

C p g(t)= f g(s)EXP s X p t ds. 
Jo 

Using (|2.1jl . one checks that C p admits a factorization 

p — OpOp, 

where 

5p = r( P + 1/2)^1. 

Consequently, it follows from Theorem 1 in JB] that the eigenvalues Ai > A2 > • • • > of C p satisfy 

as j — > 00. (2.3) 
Now the assertion follows from Theorem 1 (with ip(x) = T(p + l/2) 2 ir~ b x~ b and 6 = 2/9 + 1). □ 
An immediate consequence for fractionally integrated Brownian motions on [0, 1] defined by 

for f3 £ (0, 00) is as follows. 
Corollary 1 For every j3 S (0, 00), 

e„(y /3 )~7r-^ +1 )(/3 + lf + ^(^^) 1 /2(l g n )-(/3+^) aS n^oo. 

Proof. For p > 1/2, the Ito formula yields 

x p t={p- \) f \t - sy-lw s ds. 

2 Jo 

Consequently, 

The assertion follows from Theorem 3. □ 



Remark. The preceding corollary provides new high-resolution formulas for e n (Y 13 ) in the cases 
G (0, oo) \ N. 

One further consequence is a precise relationship between the quantization errors of Riemann- 
Liouville processes and fractional Brownian motions. The fractional Brownian motion with Hurst 
exponent p G (0, 1] is a centered pathwise continuous Gaussian process Z p = (Z£ )t e ro,i] having the 
covariance function 

EZ§Z£ = ^(s 2p + t 2p - \s-t \ 2p ). (2.5) 
Corollary 2 For every p G (0, 1), 

Proof. By we have 

1 /2 

e n (Z p ) ~ 7r-^s)(p + 1/2)" (^^) (T(2p + 1) sm^)) 1 ^^ n y P] n _ oo. 
Combining this formula with Theorem 3 yields the assertion □ 

Observe that strong equivalence e n {X p ) ~ e n (Z p ) as n — > oo is true for exactly two values of 
p G (0, 1), namely for p = 1/2 where even e„(X 1 / 2 ) = e^Z 1 / 2 ) = e n (W) and, a bit mysterious, for 
p = 0.81557. . . 

Now consider the Strassen ball U p of X p . Since the covariance operator C p satisfies C p = 
T(p + i)i? p+ i(r(p + ±)R p+ i)*, one gets 

U p = r(p + l/2)i? p+ i(B L2 (0,l)) (2.6) 

i?P+i/25 : 5 G ^ 2 ([0, 1], £ft), 1 lo5(i) 2 ^ < Hp + 1/2) 



a fractional Sobolev ball. Theorem 2 and (|2.3j) yield the solution of the entropy problem for 
fractional Sobolev balls. 



Theorem 4 For every p G (0, oo), 



■ i \ p+\ 

e n (U p ) ~ l^ 1 r(p + l/2)(logn)-^+5) 



7T 



p ^/ 2 



e n (X p ) as n —* oo. 



„ log n , 

3 Asymptotically optimal functional quantizers 

Let X be a il- valued random vector satisfying 1)1. 2j) . For every n G N, L 2 -optimal n-quantizers 
a C H exist, that is 

(E minllX-all 2 ) 1 / 2 = e n (X). 

If card (supp(Px)) > optimal n-quantizers a satisfy card(a) = n, P(X G C a (a)) > and the 
stationarity condition 

a = E(X | {X G C a (a)}), a G a 
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or what is the same 

X a = E(X | X a ) (3.1) 

for every Voronoi partition {C a (a) : a S a} (see In particular, E J" = EX. 

Now let X be centered Gaussian with dimKx = oo. The Karhunen-Loeve basis {uj : j > 1} 
consisting of normalized eigenvectors of C x is optimal for the quantization of Gaussian random 
vectors (see ^Hj). So we start with the Karhunen-Loeve expansion 

oo 

x = j2>^ /2z jUj, 

3=1 

1 /2 

where 2j =<X,uj > /A - , j > 1 are i.i.d. N(0, l)-distributed random variables. The design of 
an asymptotically optimal quantization of X is based on optimal quantizing blocks of coefficients 
of variable (n-dependent) block length. Let neN and fix temporarily m,l,ni, . . . ,n m € N with 
HJL^j < n, where m denotes the number of blocks, I the block length and rij the size of the 
quantizer for the jth block 

:= j e {l,...,m}. 

Let ay C M. 1 be an L 2 -optimal n^-quantizer for and let = Z^~> be a o^-quantization of 
Z^). Then, define a quantized version of X by 

m I 

* n ■= E E A^+^W-i)^. (3-2) 
i=i fc=i 

It is clear that card(X n (ft)) < n. Using JSU for Z^, one gets Ei" = 0. If 

^ = E Wcko,)^), 

then 

a6X™ j«j j=l fc=l 

where a = (a^\ . . . , a^) 6 XjLi a j- Observe that in general, X n is not a Voronoi quantization of 
X since it is based on the (less complicated) Voronoi partitions for Z^\j < m. (X n is a Voronoi 
quantization if I = 1 or if = ■ ■ ■ = \ji for every j.) Using again l|3.1j) for Z^> and the 

independence structure, one checks that X n satisfies a kind of stationarity equation: 

E(X | X n ) = X n . 



Lemma 1 Let n > 1. For every I > 1 and every m > 1 

m 

E||V-V"|| 2 <^A (j _ 1)menj (iV(0,/ i )) 2 + E V ( 3 - 3 ) 

j=l j>mZ+l 

Furthermore, \3.3\) stands as an equality if I = 1 (or A(j_iv + i = . . . = Aj/ for every j, I > 1). 
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Proof. The claim follows from the orthonormality of the basis {uj : j > 1}. We have 

ml i -\ 

E||X-X"|| 2 = E E Ay_ 1)I+fc E | - (ZU)) k | 2 + E a, 

J=l fc=l j>m2+l 

< E A -_ 1)J+1 E E I 4 J) - Z(i)) fe | 2 + E A, 

j=l k=l j>ml+l 

m 

= EA (i _ 1)m e n ,(^0)) 2 + E A,, 

j=l j>ml+l 



□ 



Set 

CCO^sup^/'ejfcCJVCO,/,)) 2 - ( 3 - 4 ) 
fe>i 

By (JTZJ), C(l) < oo. For every I G N, 

e nj (iV(0,/0 2 <nT 2/Z C(0 (3.5) 

Then one may replace the optimization problem which consists, for fixed n, in minimizing the right 
hand side of Lemma 1 by the following optimal allocation problem: 

m 

min{C(0 A(j-i)i+i^7 2// + A j : m, I, m, . . . , n m G N, Uf =1 rij < n}. (3.6) 

j=l j>ml+l 

Set 

m = m{n, I) := max{fc > 1 : n 1 /* \\^ 1)l+1 (n* =1 \ {j _ 1)l+1 )- l / 2k > 1}, (3.7) 

n . = n >, := [n 1 / w A;f_ 1)z+1 (n- 1 A (i _ 1) ^ +1 )-^ 2m ] J i G {1, . . . , m}, (3.8) 
where [x] denotes the integer part of x G R and 

I = l n := [(max{l,logn})' 9 ], G (0,1). (3.9) 

In the following theorem it is demonstrated that this choice is at least asymptotically optimal 
provided the eigenvalues are regularly varying. 

Theorem 5 Assume the situation of Theorem 1. Consider X n with tuning parameters defined 
in \S. 7j) - k3. y\) . Then X n is asymptotically n-optimal, i.e. 

(E ||X - X"|| 2 )V 2 ~e n pQ as n -> oo. 

Note that no block quantizer with fixed block length is asymptotically optimal (see ^J). As 
mentioned above, X n is not a Voronoi quantization of X. If a n := X n {Vt), then the Voronoi 
quantization X a " is clearly also asymptotically n-optimal. 

The key property for the proof is the following Z-asymptotics of the constants C(l) defined 
in (|3.4|) . It is interesting to consider also the smaller constants 



Q(l) := lim fc 2 /'e fc (iV(0,J0) 2 (3.10) 

fc— >oo 



(see (HZ)) ) 
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Proposition 1 The sequences (C(i))/>i and (Q(l))i>\ satisfy 

l^oo I l^oo I 1>1 I 1>1 I 



Proof. From it is known that 



liminf^ = l. (3.11) 

l— >oo I 



Furthermore, it follows immediately from (|1.7|) and (|1,8|) that 

Jim 901 = 1. (3.12) 
(The proof of the existence of lim C{l)/l we owe to S. Dereich.) For In, I € N with I > Zo ; write 

I = nln + m with n G N, m 6 {0, . . . , In — 1}. 

Since for every k £ N, 

one obtains by a block-quantizer design consisting of n blocks of length Zo and m blocks of length 
1 for quantizing N(0,li), 

e k {N{^h)f <ne [kk)/l] (N(OJ lo )) 2 + me [kl/l] (N(0,l)) 2 . (3.13) 

This implies 



k 2/i k 2/i 
C(l) < nC(?o)sup , , +mC(l)sup 



< A l ' h nC{l )+^mC{\). 



Consequently, using n/l < 1/Iq, 



and hence 
This yields 

It follows from ()3.13|) that 
Consequently 

and therefore 



C(l) A^Cjlp) 4mC(l) 
~T ~ T n + I 



C(l) A^Cjlp) 
hmsup — — < : . 

Z^oo I h 



l imsup 9M1 < liminf = 1. (3.14) 

l^oo I lo-*oo In 



Q(l) < nQ(h)+mQ{l). 



3(0 < Q(lo) ] mQ{\) 



l— >oo I In 
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This implies 

inf 9M = h (3.15) 

1 >1 l 

Since Q(l) < C(l), the proof is complete. □ 

The n-asymptotics of the number m(n, l n )l n of quantized coefficients in the Karhunen-Loeve 
expansion in the quantization X n is as follows. 

Lemma 2 (' L 12l, Lemma 4-8) Assume the situation of Theorem 1. Let m(n,l n ) be defined by \3. 7| ) 
and /TO)) . Then 

21ogn 

m{n, i n Jtn ~ — ^ — as n — > oo. 
Proof of Theorem 5. For every n G N, 

E ^-Dm^ < E + i)- 2 /< ( ^±V 

i=l J=l nj 

< 4 1 / , mn- 2 / mi (n™ 1 A (j _ 1);+1 ) 1 / m 

< 4 1 / z mA (ni _ 1) ; +1 . 

Therefore, by Lemma 1 and (|3,5|) . 



for every n € N. By Lemma 2, we have 



E||X-X"|| 2 <4 1 /'_^^A (m _ 1)m + ^ Xj 

j>ml+l 



21ogn 

mt = m(n,L n )l n ~ — - — as n-» oo. 
Consequently, using regular variation at infinity with index —b < — 1 of the function </?, 

1-6 



/2\ 

«^A ( m-i)z+i ~ mZA m j ~ f - J ^(logn) -1 



and 

_ ml^ml) 1 j 

U J as n °°' 

j>m/+l 

where, like in Theorem 1, ip(x) = l/x(p(x). Since by Proposition 1, 

hm = 1, 

n->oo l n 

one concludes ^ 

E\\X -X n \\ 2 ~ (jj ^^V(logn)" 1 as n ^ oo. 

The assertion follows from Theorem 1. □ 



Numerical and computational aspects: As soon as the Karhunen-Loeve basis (uj)j>i of 
a Gaussian process X is explicit, it is possible to compute the asymptotically optimal functional 
quantization ()3.2|) which solves the minimization problem (|3.6|) as well as its distribution and 
induced quantization error (at least for a given ??£ (0, 1)). This is possible since some optimal (or 
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at least locally optimal) vector quantizations of the N(0, /^-distribution has been already computed 
and kept off line. Let us be more specific. 

- In 1-dimension, the normal distribution N(0, 1) has only one stationary n-quantizer - hence 
optimal - since its probability density is log-concave (for this result due to Kiefer, see e.g. [3]). 
Deterministic methods to compute these optimal quantizers are based on the stationary equa- 
tion 1)3. lj) . They are very easy to implement, converge very fast with a very high accuracy. The 
Newton- Raphson algorithm is a possible choice (see |15j for details). Closed forms for the lowest 
quadratic quantization error \\Z — Z\\l2(jp> and for the distribution of the optimal n-quantization 
Z a as a function of the optimal n-quantizer a are also available in ^Sj. These three quanti- 
ties have been tabulated up to very high values of n. A file can be downloaded at the URL 
www . proba . jussieu . fr/pageper so/pages . html. 

- In higher dimension, one still relies on the stationary equation (3.1) which reads: 

E(l 0o(a) (Z)(o-Z)) =0, aea 

One must keep in mind that the left hand side of the above equation is but the gradient of the 
(squared) quantization error E||Z — Z a \\ 2 viewed as a function of the quantizer a (assumed to be of 
full size n) . A stochastic gradient descent based on this integral representation can be implemented 
easily since the normal distribution N(0, Id) can be simulated on a computer from (pseudo-)random 
numbers (e.g. by the Box-Muller method). This algorithm is known as the Competitive Learning 
Vector Quantization (or CLVQ) algorithm. It has been extensively investigated both from a 
theoretical (see e.g. ^1], pQ) and numerical (see e.g. JH] as concerns normally distributed vectors) 
viewpoints. The algorithm reads as follows: let (C(t))t>i be an i.i.d. sequence of N(0, /^-distributed 
random vectors, let (7t)t>i be a decreasing sequence of positive gain parameter satisfying J2t7t = 
+oo and J2t>i It < +°° anc ^ ^ a (0) e (K d )" denote a starting n-quantizer. Then, at time t E N, 
one update the running n-quantizer a(t — 1) := (a\(t — 1), ... , a n [t — 1)) as follows 

Competitive phase: select i(t) G argminji : \\cti(t — I) — ((t)\\ = min\\aj(t — I) — £(t)\\} 

j 

Learning phase: a i(t) (t-1) = (1 - 7*H(t) (*-!) + It C(*) 

a(t)j = oy_i(i- 1), j^i(t). 

Some further details concerning the numerical implementation of this procedure can be found in )15j . 
especially some heuristics concerning the initialization and the specification of the gain parameter 
sequence usually choosen of the form 74 = ^4^- It converges toward some local minima of the 
quantization error at a ^/77-rate. Some d-dimensional grids (d = 2 up to 10) can be downloaded at 
the above URL for many values of n in the range 2 up to 2 000. These quantizations were carried 
out to solve numerically multi-dimensional stopping time problems (pricing of American options 
on baskets, see jlfij and the references therein). 

The 1-dimensional optimal quantization of the iV(0, l)-distribution has already been used to 
produce some optimal scalar product functional quantization - i.e. based on blocks of fixed length 
1- in with some promising applications to the pricing of path-dependent European options in 
stochastic volatility models (this work is also based on results about diffusion processes from )13j). 
To be competitive with other methods (Monte Carlo, pde's) one needs to have good performances 
for not too large values of n. Within this range of values, it is more efficient to perform directly a 
numerical optimisation of (3.3) (or l)3.fi)l ) with / = 1 rather than using the theoretical asymptotically 
optimal parameters (3.7) and (3.8). 

As far as numerical implementation of functional quantization with n-varying block size is 
concerned, some first numerical experiments carried out by Benedikt Wilbertz 19^ suggest that 
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it slightly improves the scalar approach for high values of n, say n < 10 6 , simply using up to 
3-dimensional n^-quantizers with some rij not greater than 100. A similar improvement can vbe 
obtained for lower values of n (say n > 20 000) by using product quantizers made of blocks with 
mixed diemnsions (1, 2 or 3). 

Examples: The basic example (among Riemann-Liouville processes) is X 1 / 2 = W and H = 
L 2 ([0,l],dt), where 

A, = (vr(j - 1/2))- 2 , Uj (t) = V2 sin (t/Jx~) ,j>l. (3.16) 

Since for 5,p € (0, oo), 

2) 

one gets expansions of X s+P from Karhunen-Loeve expansions of X p . In particular, 

00 

X 5+ ^ = T(5 + 1) ]T ^jZjRsiuj). 

5=1 

However, the functions Rg(uj),j > 1, are not orthogonal in H so that the nonzero correlation 

between the components of (Z^ — Z(fi) prevents the previous estimates for E||X — X n || 2 given in 
Lemma 1 from working in this setting in the general case. 

However, when / = 1 (scalar product quantizers made up with blocks of fixed length 1 = 1), one 
checks that these estimates still stand as equalities since orthogonality can now be substituted by 
the independence of Zj — Zj and stationarity property (|3.1|) of the quantizations Zj,j > 1. It is 
often good enough for applications to use scalar product quantizers (see ^U], \TJ\)- If) f° r instance 
5 = 1, then 

00 

X := X 3 / 2 = ^ JXjZjRxiuj), 
3=1 

where 

= j2X~j(l - cos(t/Jx~j)). 
Note that \\Ri(uj)f = A,-(3 - 4(-1)^\/a~), j > 1. Set 



X n = J2^3Z j Rx(u j ). 

The quantization X n is non Voronoi (it is related to the Voronoi tessellation of W) and satisfies 



E||X-X"|| 2 = ^A2 J (3-4(-iy-yA J )e nj (iV(0,l)) 2 + £ A 2 (3 - ^-l)'" 1 ^). (3.17) 

j=l j>m+l 

It is possible to optimize the (scalar product) quantization error using this expression instead 
of (|3.6|) . As concerns asymptotics, if the parameters are tuned following (|3.7|) - (|3.9|) with I = 1 and 
Xj replaced by 

Vj := A 2 (3 + 4-i^Aj) ~ 3ir~ 4 j~ 4: as n — > 00, 
and using Theorem 3 gives 

(E \\X - X"|| 2 ) 1/2 ~ ( 3(12C( 4 1} + 1} \ V2 e n (X) as n ^ 00. (3.18) 
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Numerical experiments seem to confirm that C(l) = Q(l). Since Q{\) = nV3/2 (see 0, p. 124), 
the above upper bound is then 

I— — J I = 5.02357... 

4 Dimension 

Let X be a iT- valued random vector satisfying l|1.2|l . For n£N, let C n (X) be the (nonempty) set 
of all L 2 -optimal n-quantizers. Introduce the integral number 

d n (X) := min {dim span (a)) : a £ C n (X)} . (4.1) 

It represents the dimension at level n of the functional quantization problem for X. Here span(a) 
denotes the linear subspace of H spanned by a. In view of Section 3, a reasonable conjecture for 
Gaussian random vectors is d n {X) ~ 21ogn/6 in regular cases, where —b is the regularity index. 
We have at least the following lower estimate in the Gaussian case. 

Proposition 2 Assume the situation of Theorem 1. Then 

a fv\ > 1 21 °g n 

d ^ x ) ~ h i/( b -i) —g— as n^oo. 

Proof. For every n 6 N, we have 

d n (X) = mmlk>0-.e n ((g)N(0,\ j )) 2 + £ \j < e n (X) 2 1 (4.2) 

{ j=l j>k+l J 

(see HQI). Define 

c n : = mm I k > : £ Aj- < e„(X) 2 
1 i>fc+i 

Clearly, c n increases to infinity as n —* oo and by (4.2), c n < d n (X) for every n G N. Using 
Theorem 1 and the fact that ^ is regularly varying at infinity with index b — 1, we obtain 

((fc-lMcO)- 1 - £ A,~e2 n (X)~ '^A-VOogn)- 1 

j>e„+l 

and thus ^ 
Consequently, 

1 2 log n 

This yields the assertion. □ 
For Riemann-Liouville processes one concludes 



2p+l 



(see (2.3)). 



For the metric entropy problem one may introduce the numbers d n (Ux) analogously. Then, 
in the situation of Theorem 1 it is known that d n {Ux) ~ 21ogn/6 (see |12j). It remains an open 
question whether d n (X) ~ d n (Ux) ~ 21ogn/6. 
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